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JET-BOUNDARY CORRECTIONS FOR REFLECTION-PLANE MODELS IN RECTANGULAR 

WIND TUNNELS 

By BoBBBT S. SwANSON and Thomas A. Toi>i< 



SUMMAHY 

A detailed method for determining the jetrboundary corrections 
for reflection-plane models in redxmgvlar wind tunnels is 
presented. The method includes the determination of the tunnel 
span load distribviion and the derimtion of eguations for the 
corrections to the angle of attack, ihe lift and drag coefficients, 
and the pitching-, rolling-, yawing-, and hinge-moment coeffi- 
cients. The principal effects cf aerodynamic induction and 
of the boundary-induced curvature of the streamlines have been 
considered. An example is included to illustrate the method. 
Numerical values of the more important corrections for 
reflectiort^lanfi modds in 7- by 10-foot dosed wind tunnels 
are presented. 

INTRODUCTION 

The influence of the jet boundaries upon the air flow 
at and behind two-dimensional-flow models and complete 
models has been rather extensively investigated from theo- 
retical considerations. The results of several of these 
investigations are given in references 1 to 4, A few ex- 
perimental checks of the theory have been successfully 
made. The theoretical methods may be extended to deter- 
mine the influence of the jet boundaries upon the characteris- 
tics of semispan models mounted on reflection planes in 
rectangular wind timnels. One of the walls of a closed 
wind tunnel may serve as the reflection plane, as shown 
in figure 1. The jet-boimdary corrections are usually 
larger and the changes in the span load distribution are 
somewhat greater for reflection-plane models - than for 
complete models, especially with regard to the character- 
istics of the lateral-control devices. Greater care is there- 
fore required in the computations and more factors must be 
considered for reflection-plane-model corrections than for 
the usual complete-model corrections. 

The present investigation was undertaken to develop 
general methods of calculating the various corrections and 
methods of determining the changes in the span load dis- 
tribution caused by the jet boundaries. Nimierical values 
of the more important corrections were calculated for a 
series of representative models mounted in 7- by 10-foot 
closed rectangula;r wind tunnels. The numerical values are 
presented in the form of graphs and empirical equations in 
a separate section of the report, in order that the values may 
be obtained without referring to the detailed calculation 
procedure. Tables of the numerical values of the jet- 
boundary-induced upwaah velocity for 7- by 10-foot closed 



wind tunnels are included and should be used if it is desired 
to compute the corrections for models having unusual pro- 
portions. The complete calculation procedure is illustrated 
in detail by an example. 

The basic method used to determine the jet-boundaiy 
corrections is to determine the increments of aerodynamic 
forces and moments acting on a model which is twisted by 
the amount of the boundary-induced upwash angle. Methods 
of calculating the boundary-induced upwash angle along the 
model span and chord and methods of calculating the various 
jet-boimdary corrections, accounting for the principal 
effects of aerodynamic induction, are presented in separate 
sections in the present report. 

The formulas and corrections presented apply to com- 
plete models for which the spans are twice the spans of the 
reflection-plane models. If a model of only the outer wing 
panel is tested, the measured characteristics will be for a 
model of the aspect ratio, taper ratio, and lateral-control- 
device span ratio actually tested. Additional plan-form 
corrections — that is, the usual aspect-ratio and taper-ratio 
corrections plus correctioms for the ratio of the lateral- 
control-device span to the wing span, reference S^must 
therefore be made to determine free-air data for the actual 
airplane from the corrected data for the model. 

SYMBOLS 

r I 'circulation strength of vortex 

Ci measured lift coefficient 

ACi; increment of lift coefficient 

^^*».e. correction to lift coefficient due to streamline 

curvature 
Ct section lift coefficient 

Cn section normal-force coefficient 

Cft section hinge-moment coefficient 

AL, increment of lift at any section 
(7, measured roUing-inoment coefficient 

Che corrected rolling-moment coefficient 

A(7i increment of rolling-moment-coefficient correction 
due to jet boundaries other than reflection 
plane, based on a reflection-plane rolling- 

2Aa^ 

moment coefficient of 

Oijk rolling derivative due to deflection of lateral- 
control device (referiance 5) 
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half of increment of rolling-moment-coefficient 
correction due to reflection plane based on unit 
free-air rolling-moment coefficient , 

aerodynamic-induction factor used in determining 
AC, 

angle of attack 

correction to angle of attack 

increment of iaduced drag at any section 

correction to induced-drag coefficient 

correction to induced-yawing-moment coefficient 

[(AC7.,)x+(AC7„),+(Aa,),4-(Aa.,),+ (Aa.,)J 

increment of yavring-moment-coefficaent correc- 
tion due to reflection plane 

increment of yawing-moment-coefficient correc- 
tion due to boTmdary-induced aileron upwash 
and Tsdng loading 

increment of yawing-moment-cbeffident correc- 
tion due to boundary-induced Tving upwash and 
tunnel aileron loading 

increment of yawing-moment-coefficient correc- 
tion due to boundary-induced aileron upwash 
and flap loading 

increment of yawing-moment-coefficient correc- 
tion due to boundary-induced flap upwash and 
tunnel aileron loading 

correction to pitdung-moment coefficient due to 
streamline curvature 

correction to hinge-moment coefficient 

correction to hinge-moment coefficient at any 
section 

increment of hinge-moment correction at any 

section 
air density 

free-stream velocity, paxallel to X-axis 

dynamic pressiure Qp^*^ 

induced vertical velocity, parallel to Z-axis 
distance parallel to X-axis 
distance parallel to Z-axis 
distance parallel to r"-axis 
centroid of spanwise load 
spanwise position of trailing vortices 
effective height of wing above tunnel center line 
total wing span (twice span of reflection-plane 
model) 

span of aileron on semispaa model 

span of flap on semispah model 
total wing area (twice area of reflection-plane 
model) 

area of aileron on semispan model 
area of flap on semispan model 

aspect ratio 

taper ratio, fictitious chord at tip divided by 

chord at root 
chord at any section ' 
mean chord 



c, chord at plane of symmetry 

r radiiis of curvature of streamlines 

h tunnel height 

a tunnel breadth 

00 slope per radian of section lift curve 

01 slope per radian of lift curve of finite-span wing 
Pmax maximum ordinate of jet-boundary-induced ellip- 
tical load 

F hinge-moment correction factor for jet-boundary- 

induced elliptical load 



Subscripts: 


w 


wing 


f 


flap 


a 


aileron 


b 


overhang balance 


av 


average 


t 


tunnel 


total 


total 


c 


corrected 


r 


reflection 


p.p. 


principal part 


8. p. 


supplementary part 


On. 


for pitching moments 


e 


effective 


s.e. 


streamline curvature 


l.l. 


lifting line 




spanwise location of trailing vortex 



The axes used are defined in flgm-e 1. All loading and 
botrndary-induced upwash-velocity parameters with primes 
are based on lift or rolling-moment coefficients not equal 
to unity. 

BOUNDABY-INDUCED UPWASH 

THEORY 

General problem. — The general problem to be solved in 
determining the jet-boundary corrections for a complete 
model in a wind tiumel is the determination of the total up- 
wash velocity induced by the jet boundaries. The special 
problem for a semispan model mounted as a reflection-plane 
model to simulate a symmetrically loaded complete model is 
the determination of the total boundary-induced upwash 
velocity minus the induced upwash velocity due to the re- 
flection of the semispan model. The problem of determining 
the boimdary-induced upwash velocity duo to imsymmetrical 
loading devices, such as lateral-control devices .on reflection- 
plane models, is one of determining not only the total 
boundary-induced upwash velocity, as for complete models, 
but also an additional correction due to the nonexistence of 
the reflection wing. 

Use of images.— ^The known boundary conditions to be 
satisfled are zero normal velocity for closed wind tunnels and 
constant pressiire for open wind tunnels. The boundary 
conditions for a closed rectangular tunnel may be satisfied by 
a doubly infinite system of images (reference 1). Figure 1 
shows the three-dimensional image arrangement that satis- 
fies the boundary conditions for a semispan model mounted 
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in a closed recfcangiilar wind tunnel. The model is mounted 
on the XZ-plane — or left wall, looking downstream — and 
located in the Xy-plane. The reflection wing is shown ia 
phantom and lies along the negative y-axis. It may be noted 
that this image arrangement is the same as that for a com- 
plete model of the same semispan ia a tunnel of the same 
height and twice the width. The images of the wing are 
represented in this figure as simple horseshoe vortices of 
semispan \yi\. Any actual span load distribution may be 
constructed to any desired degree of accuracy from a com- 
bination of several horseshoe vortices. The three-dimen- 
sional image arrangement is necessary only when it is de- 
sired to compute the boundary-induced upwash velocity 
behind tlio lifting line, the streamline curvature, or the 
boundary-induced upwash velocity for models with excessive 
sweepback. 

The boundary-induced upwash velocity at the lifting line 
may be determined from a two-dimensional image arrange- 
ment satisfying the boundary conditions at infinity as shown 
by Prandtl. Figiu-e 2 shows the two-dimensional image 
arrangement satisfying the boundary conditions for a single 
counterclockwise trailing vortex and its reflection (clock- 
wise) located at a distance d above the tunnel center line 
and at distances equal to yi and —yi from the reflection wall. 
The single trailing vortex and its reflection represent the 
trailing vortices of a simple horseshoe vortex with semispan 
equal to \yi\. 

CALCITLATION METHODS 

Preliminary calculations.— The calculations of general 
curves of boundary-induced upwash velocity for varioiis 



imago arrangements (figs. 1 and 2) of simple horseshoe 
or trailing vortices will considerably simplify the labor 
involved in determining the boundary-induced upwash 
velocity for any given model. The boundary-induced 
upwash velocity behind the liftiag line for two values of 
vortex semispan (image arrangements of %. 1) was calcu- 
lated by the methods described in reference 2. The results 
of the calculations are presented in figure 3 and table I. 
These calculations apply either to a reflection-plane model 
in a 7- by 10-foot closed wind tunnel or to a complete 
(symmetrically loaded) model in a 7- by 20-foot closed 
wind tunnel. 

The boimdary-induced upwash velocities at the lifting 
line were obtained by computing the combined effect of 
enough of the images, corresponding to the arrangement 
of figiu:e 2, to give values accurate to the foiu^ decimal 
place. The results are given in figiu-e 4 and table II. 
These values apply to complete models mounted in 7- by 
20-foot closed wind tunnels, as weU as to reflection-plane 
models mounted in 7- by 10-foot closed wind txumels. 

The method used to determine the boundary-induced 
upwash velocity at the lifting line for any given image 
arrangement is to break up that image arrangement into 
certain groups, usually vertical rows of images, for which 
simple summation formulas are available. The sum of the 
effects of each of these groups may then be determined. 
The summation formulas for vertical rows of vortices 
extending from th'e FZ-plane to infinity in one direction 
were developed in references 3 and 4. 
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FiODBE I.— Tbrce-dlmesslonal anangemeDt of the dooblr infinite Image pattern satldTlng the bonndary conditions for a reilection-plane model of a wing on tbe center line of a dosed 

rectangular wind tunnel. 
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Tlpwasli velocity for nonaniform span loading. — ^Upwash 
velocity for any nonuiuform span loading may be approxi- 
mated to any desired degree of accuracy by breaking down 
the actual loading into several steps over each of which the 
loading is assumed uniform. The boundary-induced upwash 
velocity may then be determined as the sum of all the com- 
ponents of upwash velocity due to all the stepwise incre- 
ments. Numerical values for the upwash velocity may be 
taken directly from the tables rather than from the figures, 
provided that increments are taken at J4-foot values of yi. 
If the tunnel walls appreciably alter the span-load distribu- 
tion, as they usually do for a model with a lateral-control 
device having a relatively large span, the actual tunnel span 
load distribution should be used instead of the theoretical 
free-air span load distribution. Methods of approximating 
the tunnel span load distribution wiU be presented later in 
this report. A three- or four-step appro:dmation to the tun- 
nel span-load curve is usually necessary for asymmetrical 
load conditions. Calculations indicate that very lai^e errors 
are introduced by using single uniform loadings for the 
asymmetrical conditions. For synunetrical load conditions, 
however, a single-step approximation is usually satisfactory 
if the span over which the imiform load is assumed to act 
(called effective span) is properly chosen. 

The span-load parameter that will be used in the compu- 
tations is cciAIICl, which is eqmvalent to 2TA/bVCz, (and 
to La as used in reference 6). The upwash angle in radians 
for unit lift coefficient is w/VOz. for small angles. The 
formula for determining the local upwash' angle is then 



(1) 



where 0)^^ is obtained from figure 4 or table Tl and 

A(eciA/bCL)vi is proportional to an increment of load ex- 
tending from the reflection plane to yi. In other words, 
A(cctA/bOi,)f^ is proportional to the strength of the trailing 
vortex assimied to leave the wing at yi. 

Boundary-induced upwash angles are given in figure 5 for 
unit lift coefficient for a 7-foot semiapan model of aspect 
ratio 6 and taper ratio 0.5 and for a unit flap lift coefficient 
for two ratios of flap span to wing span (called flap-span 
ratio). The actual span loading is represented by a seven- 
step approximation. It may be seen from figure 6 that, if 
the proper value of the effective semispan is used, the up- 
wash angle may be determined satisfactorily by the use of a 
simple uniform load. The effective span is, of course, de- 
pendent upon the particular model-tunnel configuration. 
Computations for several representative reflection-plane 
models in 7- by 10-foot closed wind tunnels indicate that 
good accuracy in the determination of the boundary-induced 
upwash angle is obtained with the following- ratios of ef- 
fective span to actual span, bjb or bfjb, provided the model 
semispan is between 6 and 8 feet (the usual values): 

Wing: 

Rectangular 0. 03 

Taper ratio, about 0.5 i 0. 88 

Taper ratio, about 0.25 - 0. 83 

Partial-span flap: 

bfib greater than 0.6 — 1. 00 

b/fb less than 0.6 , 1. 30 



3 



G 



3 



3 



G 



■I- 



3 



■f 

I 

1- 



3 
3 



4- 



3 I G 

h 

-y-r 



-a > 



G 



3 



3 ^ {4- G 



U— y.- 



I t 



3 I G 



3 



^ 

3 j G 

I 



G 



I 

•4- 



3 1^ G 



3 

3 



G 



r 



3 



G 



3 



G 



G 



FiQxmE 2.— TTwo-dlmeiuioiiBl airanEomeiit of the infinite Image p&ttom satisfying the bbnndary conditions for & single trailing vortex and its refleotlon located above the center Uno of a 

closed rectangular wind ttmneL 
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O 2 4 '6 a - 10 



Disfonce from pkmo of symmetry, y, ft 

FlOOTE t-Bonndary-lndnoea upwash vetodty at the mUng Una due to a stngle ooonterclockwlse trafllng 701161 looated at various distances ft from the reflection plane 

dosed leotangnlar \rind tunnels. 
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The span-load parameter that ■will be used in the com- 
putations for the asymmetrical loading condition is cciA/bCi 
(or 2VAJbVCi). The upwash angle, however, is most con- 
veniently expressed in terms of a parameter that represents 
the reflection-plane loading for a rolling-moment coefficient 
greater than unity, as foUovre: 

[y(C7,,-r2A(7o] ^ (3\^ [6(Ci.+2A(7,,)] 

The primes indicate that the expressions are not true param- 
eters, because the values are not for unit rolling-moment 
coefficient but for a rolling-moment coefficient equal 

to 1+-^'- 

CORRECTIONS 

GENERAL 

The jet-boundary corrections may be divided into two 
groups. The first group consists of all corrections to be 
applied to a symmetrically loaded model; that is, all forces, 
moments, and air-flow conditions acting on the reflection- 
plane model are reflected identically with respect to reflec- 
tion plane and thus the measured model characteristics axe 
for a symmetrical model. The measured lift, drag, and 
pitching moment of the model are thus exactly one-half 
those for a complete model mounted in a wind tunnel of 
the same height and twice the breadth as the original tunnel, 
and the boimdary-induced upwash velocity is the same as 
for the complete model in the larger tunnel. The second 
group of corrections is for the asymmetrically loaded con- 
dition, such as the loading due to the deflection of a lateral- 
control device. The loading due to the lateral-control de-. 
vice is reflected into the reflection plane just as it was for the 
symmetrical case; but the reflection is undesirable in this 
case and must be corrected for, as it would not be present 
on a complete model. Also, the absence of the other wing 
(the reflection wing) causes the measured rolling and yawing 
momenta to be too lai^e, because the load due to aerodynamic 
induction existing on the other wing of a complete modd 
will be absent from the measured values of a reflection- 
plane model. 

The corrections will be determined with the proper sign in 
order that they may be added to the measured values for 
a closed-throat wind tunnel. 

STAIMETBICAL LOADING CONDmONS 

Span load distribution. — The correction to the span load 
distribution need not be determined imless stalling tests or 
actual span-load-distribution tests are made. Calculations 
for a few reflection-plane models of usual size in 7- by 10-foot 
closed wind tmmels indicate that the wing span load dis- 
tribution is altered by the tunnel walls by an amount 
equivalent to a change in taper ratio of about 0.05; that is, if 
the geometrical taper ratio of the model is 0.50, the wing 
span load distributionTn the tunnel corresponds to a wing 
having a taper ratio of about 0.45. The changes in flap 
span loading axe somewhat greater than the changes in 
wing span loading. The usual effect of the tunnel walls on 
the flap span loading is to increase the relative loading 



over the unflapped portion of the wing and to reduce the 
relative loading over the flapped portion of the wing. It 
should be remembered, however, that the type of change in 
span load distribution caused by the jet boimdaries is 
entirely dependent upon the model-timnel configuration and 
that other model-tunnel arrangements might produce effects 
opposite to those just indicated for a reflection-plane model in 
a 7- by 10-foot closed wind tuimel. 
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FiGDBE 6.— Values of ejllb to be tucd In conrertliig tho load parameter cei/c.a of rcfcrenco 6 
to tbe-Joad parameter cetAlba. 

The span load distribution of the wing in free air must bo 
determined to evaluate the change in loading due to tho 
tunnel walls. The free-air span loading for symmetrical load 
conditions may easily be obtained from the tables of refer- 
ences 6 and 7 for several different wing-flap combinations. 
For other flap arrangements or any initial ^ving twist, tho 
influence lines of reffflrence 6 may be used to estimate the 
shape of the load curves. The actual load curve may be 
determined from the condition that the area imder the curve 

of cciA/bGz, plotted against ^ is equal to unity or, if 

plotted against y, is nmnerically equal to 6/2; that is, tho 
average ordinate is equal to- unity. 

The increment of boundary-induced load corresponding 
to a tunnel lift coefficient of unity is obtained by asstuning 
that the wing is twisted an amount equal to the boundary- 
induced upwash angle w/VOi.. The boimdary-induced load is 
calculated from the influence lines of reference 5. The in- 
fluence lines give the load at a particular spanwise station 
for imit changes in angle of attack extending various dis- 
tances inward from the wing tip. Values of c,A/b for the 
wings of reference 5 are given in figure 6 in order that the load 
parameter cci/c,a used in reference 6 may be converted to 
cciAP)oc. 
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The incroment of boundary-induced load determined in 
this fashion may be added to the free-air span-load curve to 
obtain a first approxiination to the tunnel-load curve cor- 

responding to a tunnel lift coefficient of 1+-^} where 

^OlICl is equal to the average ordinate of this increment of 
boundary-induced load. Because the tunnel lift coefficient 

1 is greater than unity, a second approximation to the 

incroment of boundary-induced load based on a lift coeffident 

of 1 and the new load distribution must be determined, 

and then a third approximation must be made, and so forth. 
In order to avoid the necessity for using successive approxi- 
mations in this fashion, it may be assimied that th« values of 
w/VCz, used for the first approximation need only be multi- 
plied by constants determined from the increase in lift co- 
efficient for each of the remaining approximations; that is, 
it may be assumed tliat the change in the shape of the span- 
load curve would not change the shape of the w/VOi. curves. 
It was shown that the values of- the w/VOi, curve may be 
computed with satisfactory accuracy when a uniform loading 
over an efi'ective span is assumed and, Inasmuch as the shape 
of the tunnel-load curve does not change appreciably, this 
assumption of unchanged wIVGl curves is reasonable. The 
increment of boundary-induced load corresponding to the 
approximation based on a free-air unit lift coefficient (if 
the shape of the wIVOl curve is unchanged) is simply equal 



to the values for the first approximation times 



" Or. 



In 



order to obtain the loading in the tunnel corresponding to a 



lift coefficient of 



A^ 



'^Q'> add to the free-air load 



this 



approximation. Divide all these values by 1-f- 



AOi. 

Gr. 



AGi. 

to obtain the tunnel-load curve for unit lift coefficient. 

This method of estimating the tunnel span load distri- 
bution takes into account the main effects pi aerodynamic 
induction. The method is not exact because the jet- 
boundar^-induced upwash angle w/VOl is calculated approx- 
imately. If desired, the upwash angle corresponding to the 
tunnel-load carve previously determined can be obtained 
with great accuracy by using many steps in the stepwise 
distribution to represent tiie tunnel-load curve. The 
calculations can then be repeated with the new values of 
boundary-induced upwash angle. The process could be 
repeated xmtil the exact tunnel loading, insofar as lifting- 
line theory applies, is obtained. It seems, however, that 
the process is so rapidly convergent that the span loading 
calculated from the approximate upwash angles is usually 
satisfactory. 

As a check on the convergence, the tmmel span load 
distribution for a large-span elliptical wing in a circular wind 
tmmel was computed by the method previously described and 
the result was compared in figure 7 with the more exact 
calciilation of the tunnel span load distribution made by 
MiUikan (reference 8). It may be noted from figure 7 that 
the effect of the tunnel walls is opponte in this case to the 
effect already described for models of usual size in 7- by 
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FiouBE 7.— Span load dlstribntlons tor an elUpUcal wing in free air and in a closed olrcular tnnneL Wing span is 0.9 times the tunnel diameter. 
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20-foot closed wind tuimels or reflection-plane models in 
7- by 10-foot closed wind tunnels. Influence liaes, sinular 
to those given in reference 5, were determined for an elliptical 
wing for use in calculating the boundary-induced load 
increment. 

A comparison of the final timnel-load carve with the 
original (free-air) load curve indicates the change in span 
loading caused by the tunnel walls. The increment of 
load due to streamline curvature may also be added to the 
original load curve. Because the original curve was ob- 
tained, however, from the lifting-line theory rather than 
from the lifting-surface theory, such an additional step 
would seem an undue refinement. 

Chordwise load distribution. — ^The chordwise load and the 
chordwise load distribution at each section are altered 
by the jot boundaries. The main portion of this change 
in load is corrected for by the usual induced angle-of-attack 
correction for the upwash angle at the lifting line. The 
curvature of the streamlines caused by the jet boimdaries 
effectively changes the airfoil camber, which results in a 
further change in the chordwise load (and the chordwise 
load distribution). The corrections due to the change in 
effective camber may be applied partly as an increased 
angle-of-attack correction and partly as a correction to the 
Jif fc, the pit<;hing moment, and the hinge moment. 

The general characteristics of the increment of load 
due to boundary-induced streamline curvature may be esti- 
mated from thin-airfoil theory. The shape of the boundary- 
induced streamlines is, to a first approximation, circular 
because the boundary-induced upwash angle varies al- 
most linearly along the chord unless very wide-chord models 
are used (fig. 3). 

The chordwise load for an airfoil with circular camber 
may be broken into two components. One component 
corresponds to a loading of the flat-plate type, which is 
similar to the loading due to an angle-of-attack change 
(also called additional-type loading). The magnitude of 
the load is determined from the product of the dope of the 
lift curve and the boundary-induced increase in the angle of 
inclination of the tangent at the 0.50c point because, for 
circular camber, the curve at this point is parallel to the 
chord line connecting the ends of the mean line. Inasmuch 
as this component of load is similar to the load resulting 
from a simple angle-of-attack change, it may be applied 
as an additional angle-of-attack correction. The other 
component of load is eUiptically shaped and its magnitude 
is determined by the product of the slope of the lift curve 
and the angular difference between the zero-lift line and the 
chord line or the 0.50c-point tangent line. The zero-lifL 
line, for an airfoil with circular camber, is determined by the 
angle of inclination of the tangent of the 0.75c point. The 
lift, the pitching-moment, and the hinge-moment corrections 
are a result of this elliptical component of load. 

The location of the lifting line for a plain airfoil may be 
assumed to be at the 0.25c point and the boundary-induced 
upwash angle is computed by assuming that the total lift is 
concentrated at the lifting line. The location of the lifting 



line for a flapped wing wiU lie somewhere behind the 0.25c 
point, depending upon the flap characteristics. The location 
of the lifting line determines the magnitude and direction of 
the flat-plate type of load. The two components (flat-plate- 
load and elliptical load) are equal and positive if the lifting 
line is located at the 0.25c point. Each component may be 
expressed as 

Ac,=0.25 ^«0.25aocCi 

^ ox 

where r is the radius of curvature of the streamlines. If the 

lifting line is at the 0.50c point, the flat-plate component of 
load is zero. The elliptical component is positive ancl equal 
to 0.25aoc/r, because it is independent of the location of the 
lifting line. If an extensible flap is used, the magnitude and 
the distribution of the chordwise load must be calculated as 
though the chord of the airfoil were increased. Because the 
results of the tests are usually based on the original chord, 
the final correction must also be reduced to a coofficiont 
based on the original chord. 

No correction wiU be appKed directly to the chord^v^so 
load distribution but the angle of attack, the lift, the pit<?hing 
moment, and the hinge moment will be corrected to accoimt 
for the altered load distribution. 

Angle of attack. — ^The main portion of the angle-of-attack 
correction is due to the jet-bomdary-induced upwash angle 
at the Hftiog line. The problem of finding the angle-of- 
attack correction is, basically, the determination of the angle 
of attack that^ the model woiUd require in free air to have 
the same lift as in the wind tunnel. The correction angle is, 
then, the difference between the free-air angle of attack and 
the tuimel angle of attack. If the tunnel span load distribu- 
tion is determined, the angle-of-attack correction duo to the 
boimdaiy-induced upwash at the lifting line is given as 



. AO; 57.3^ 



(3) 



where ACJGl is the increment of boundary-induced load for 
a tunnel lift coefficient of unity as determined for the span- 
load calculations. 

If the tuimel span load distribution is not determined, the 
angle-of-attack correction may be calculated by an alternate 
method that gives values almost identical with those of 
the method just described. For the alternate method 
the boundary-induced upwash angle is weighted according 
to the wing diord at each section and is then averaged 
across the span. The formula is 



. 57.30; w 



(4) 



The increment of additional load caused by streamline 
cmvature is dependent upon the relative distance between 
the lifting line and the 0.60c point as indicated in the section 
on chordwise load distribution. In the case of the Aving, the 
lifting line may be assumed to be located at the 0.25c point ; 
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thus 0.50— ^'=0.25. Becaiise the lifting line due to de- 
flection of a partial-span flap is usually located very near 
the 0,60c point, it generally is not necessary to apply a 
correction to the angle of attack for this case. The general 
equation for the correction to the angle of attack is 



^(w)2A 
where the quantity — -j- is equivalent to 



c'dy (5) 



Values of 



da; 



may be obtained from table I for 



2/1 = 3 feet and for 7/1=6 feet. The correction angle determined 
by equation (5) should be added to the lifting-liue correction 
angle as computed from equation (3) or from equation (4). 

lift. — ^The measured lift in a closed wind tunnel is greater 
than it would be in free air even though the complete angle- 
of-attack correction is applied. The increase in lift is due 
to the elliptically shaped increment of chordwise load 
caused by streamline curvature. This elliptically shaped 
increment is applied as a correction to the lift rather than as 
a correction to the angle of attack in order to correct the 
maximum lift coefficient. This increment is determined 
from the slope of the lift curve and the difl'erence in jet- 
boundary-induced upwash angle at the 0.50c point and at 
the 0.76c point. The integral 



da; 



^Y<^dy calculated for the Aa,.e. correc- 



tion may nlso bo used for the lift correction 



ACi =- 



46 ^0 bx h 



<?dy 



(6) 



where Oi is used instead of Oq to account approximately 
for aerodynamic induction. 

If the lifting line can be assumed to be located at the 
0.25c point, as in the case of the wing, equation (6) may be 
written for the wing lift correction as 



(A<:?£..,.)„=-(A«..c.)„ 



57.3 



(7) 



A small correction to the lateral center of pressure, which 
is detormiaed from the rolling moment, may be obtained by 
performing a moment integration of the streamline-cxuvature 
load, as 



A(7,. 



0 



'(wJ: 



to 

da; h 



It should be noted that this increment is based on the 
complete span 6 rather than on the model span 6/2, 

Pitching moment, — The pitching-moment coefficient must 
be corrected for the elliptical component of the boundary- 
induced load. The correction is 

da; 



AC„=- 



46c„ 



(?dy (8) 



where the factor ""^'''^J accounts approximately 

for the effects of aerodynamic induction and reduces to unity 
if the moments are taken about a line through c»/4. The 
distance is measured between the line about which the 
moments are computed and the midchord at the spanwise 
centroid of the elliptical load y,,^, where 



2A 
b 



bx 



Jo 6 oa; 



(9) 



Downwash behind wing. — The correction to the dowawosh 
beMnd a reflection-plane model is determined from the 
boimdary-induced upwash curves previously computed. 
The general methods of reference 2 should be used to deter- 
mine the downwash-angle correction. 

Drag. — The induced-drag correction is determined from 
the generalized KuttarJoukowsM law. The boundary- 
induced upwash angle and the tunnel span loading ore used 
in the computation of the induced-drag correction. The up- 
wash angle at each section must be multipHed by the loading 
at that section and the result integrated mechanically. In 
order to establish the method and to determine the propor- 
tionality constants, the integration formula wiU be 
developed. 

The increment of induced drag at any section due to the 
tunnel walls is 

AD,=/ywT dy (10) 
and the increment of induced-drag coefficient for the wing is 



AG, 



, 4_ f W2 



wF dy 



(11) 



where the product wT is the sum of all component products 
of wing and flap upwash velocity and circulation 
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Flap- or oileron-chord ratio, '^O'"^ 

FicvRB 8.— Variation ot the binge-moment conectton lactor Ibr dUptical load F Trltb the flap- 
or alleron-cbord ratio c//c or cjc and the balance-chord ratio eile/or eijc 

A correction to the lateral center of pressure of the induced- 
drag force may be obtained by integrating the increment of 
drag given, in equation (10) for the yawing-moment 
coefficient 



AG, 



2 r*/2 



wTydy 



(13) 



Hinge moments. — ^The measured hinge moments of the 
high-lift device and lateral-control device should be cor- 



rected for the component of eUiptical load caused by the 
curvature of the streamlines. The hinge-moment correction 
is determined from an integration of the moment about the 
hinge axis of this load on the high-lift device or on the con- 
trol smrface. The integration must be performed over the 
entire surface — both chordwise and spanwise — ^because, in 
general, the correction varies along the span. The chord- 
wise integration can be performed analytically because the 
shape of the load is known to be nearly elliptical. The 
increment of load at any section is determined from the 
expression for the area of an ellipse as 



(14) 



where Pmax is the maximum ordinate of the elliptical load. 
The increment of load is also equivalent to 



8A 



2A, 



(15) 



The increment of hinge moment AH, at any section is 
obtained by a moment integration about the hinge axis of 
the part of the elliptical load over the movable surface, as 

Aff,=fp_J7l-(^J(^+^-l)d(^^) (16) 
ax-y/l— is the ordinate of the elliptical load. 



where P, 



The value of Pmu is obtained by solving equations (14) and 
(15) and is substituted in equation (16) to give 



AH,=^ 



'b 



c/2 



0<i) 

where the integral 

wiU be known as the hinge-moment correction factor for 
eUiptical load F and the integration is carried across the 
flap and across any overhang type of balance. The results 
of the integration are presented in figine 8 as a fimction of 
the ratio of the flap chord to the airfoil chord C//c or of the 
aileron chord to the airfoil chord cJc, and of the overhang 
balance chord ratios Cj/c/ and Cj/ca. The corrections to the 
flap hinge-moment coefficient may now be expressed as 



AG, 



GzflJ> 
'SvASfCf 



dy 



(17) 



where the spanwise integration across the flap must be 
performed mechanically. The correction to the aileron 
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binge-moment coefficient may be determined by performing 
the integration over the limits of the aileron span rather 
than of the flap span, as indicated in equation (17), and by 
using Sa and Ca. It should be noted that is a fimction of 
only Cf/c (or Ca/c) and c^/c/ (or cjca) and will therefore have 
the same value at all sections of constant-percentage- 
chord flaps or ailerons. The efiPect of aerodjoiamic induc- 
tion on the hinge-moment correction due to streamline 
curvature is small and will be neglected. 

ASYMMETRICAI. LOADING CONDITIONS 

Span load distribution. — The jet boimdaries have a pro- 
nounced effect upon the span load distribution of asym- 
metrical load devices on reflection-plane models. In order 
to determine the rolling- and yawing-moment corrections,, 
part of the computations to determine the tunnel span load 
distribution must be made. The actual distribution may 
be obtained by a small amount of additional work and a 
more accurate estimate of the rolling-, yawing-, and hinge- 
moment corrections is then possible. 

The tunnel span load distribution is determined by adding 
to the free-air load the increment of load due to the reflection 
plane and the increment of load due to the other jet boun- 
daries and by then reducing this total load to that for a 
rolling-moment coefficient of imity. The influence lines of 
reference 5 may be used to estimate the free-air load, and 
the reflection-plane load increment. It should be noted 
that the reflection-plane load increment is simply the load 
induced on the reflection wing for an asymmetrical load on 
the real wing; that is, the load curve for a reflection-plane 
model in free air (no jet boundaries- except the reflection 
plane) is obtained by adding the free-air load at —y to that 
aty. 

The load parameter ecilc^ of reference 5 should be 
changed to cCiA/bOi^ for convenient use in the computation 
of jet-boundary corrections. The conversion may be made 
as 

CC,A_ CCi l^lc ^ jgj 



60/; C,a b I a. 

where values of c,A/b are given in figure 6 for the wings of 
reference 5 and values of Ctja may be determined from 

figure 16 of reference 5 — that is, ~=0.5-j^- The conver- 
sion may be made graphically from the condition that the 
moment of the area under the curve of cciA/bOi^ against 

^ is equal to 4.0 or, if plotted against y, is numerically equal 

to 4.0(6/2)' or b\ 

The increment of load due to the jet boundaries (other 
than the reflection plane) is obtained from the influence 
lines of reference 5 by the same general methods used in 
determining the increment of boimdaty-induced load for the 
symmetrical loading condition; that is, the -wing is assumed 
to be twisted by the amount of the boundary-induced 
upwash "angle and the corresponding increment of load is 
obtained from the influence lines. 



The boimdary-induced upwash angle should, strictly 
speaking, be determined by a process of successive approxi- 
mations because it also depends upon the shape of the 
tunnel span-load curve. It is usually satisfactory, however, 
to represent the reflection-plane loading by a three- or fmu-- 
step approximation, and to calculate the corresponding 
boimdary-induced upwash angle by the methods suggested 
in the section on boundary-induced upwash velocity. The 
increment of load calculated from this boundary-induced 
upwash angle will correspond to a rolling-moment coefficient 



of 1-f 



2A0i 



-) which is the reflection-plane rolling-moment 



coefficient for imit free-air rolling-moment coefficient. (The 
use of these factors -wUl be made clearer in the section on 
rolling-moment corrections and in the iQustrative example.) 
The increment of load must therefore be increased to 
correspond to the rolling-moment coefficient occnrring in 
the tunnel for -unit free-air rolling-moment coefficient 
((7|„=1.0). It will be sho-wn later in the section on rolling- 
moment corrections that the tunnel rolling-moment coeffi- 
cient for unit free-air rolling-moment coefficient is equal to 

, 2A0,, 



1 + 



2A^ 



where AC, is the moment of the increment of boundary- 
induced load corresponding to the reflection-plane load. 
This increment must therefore be multiplied by 



2AG, 



1+ 



before it is added to the reflection-plane load. The tunnel 
span load distribution for unit tunnel rolling-moment 
coefficient is obtained by so reducing the ordinates of this 
curve that the moment is equal to 4.0(6/2)* or (the same 
thing) by multiplying by the rolling-moment-coefficient 
correction 

AO, 



1^ Gt. 



The timnel span load distribution, as weU as the roUing- 
moment-coefficient correction, has thus been determined. 
The explanation of the determination of the rolling-moment 
correction will be ^ven in the next section in some detail 
to explain further the method of determining the tunnel 
loading and to present alternate methods of determining the 
rolling-moment correction that do not require the determi- 
nation of the tunnel loading. 
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Aileron span ratio, ■^Jg 



2ACi; 



FiQxmE 9.^Va]tiee of the &ctor iH — used to accoimt Ibr the effect of the lefleotlon plane on the ndlln&'moment 
coefficient for idleiona extending inward from the tips of the wings of reference 6. 



Eolling moment. — The correction to the rolling moment 
will be determined in two parts. The first part of the 
correction is caused by the absence of the reflection wing 
and the load increment due to the reflection plane. This 
part of the correction depends not upon the model-tunnel 
configm*ation but only upon the characteristics of the model 
itself; consequently, it was possible to calculate the cor- 
rection increment for several wing-aileron combinations 
from the data of reference 6. The aileron span load distri- 
butions were mechanically integrated to determine the 



moment of the load on the absent wing AOtJCi^ in terms of 
the free-air moment of the total load. Not only because the 
reflection wing is nonexistent but also because an equal load 
is induced on the real wing by the reflection image, t\vice 
this absent-wing moment must be applied as a correction. 
The correction is presented in figure 9 in an easily used form, 



as a fimction of the aileron span ratio — that is, 



ratio of aileron span to wing semispan- 
ing inward from the wing tips. 



-for ailerons extend- 



/ 
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The second part of the roUing-momeiit correction. A<7, re- 
sults from the moment about the plane of symmetry of the 
boundary-induced load. The method of calculating the 
boundary-induced load has already been explained in con- 
nection with the determination of the tunnel span load 
distribution. If the tunnel span loading is not determined, 
this second increment may be c^culated from simple strip 
theory — neglecting aerodynamic induction, that is, the 
effects of the velocities induced by the trailing-vortex 
system — and multipKed by a factor A/(A+J) to account 
approximately for the effects of aerodynamic induction. 

The value of J depends upon the distribution of boundary- 
induced upwash angle and the model taper ratio but is 
practically independent of aspect ratio and only slightly 
dependent upon the slope of the section lift curve. Figure 10 
gives some values of J as a function of taper ratio for three 
distributions of boundary-induced upwash angle. The 
corresponding upwash-angle distribution is also given in 
figure 10 and corresponds to various aUeron-span ratios for 
reflection-plane models in 7- by 10-foot closed wind tunnels. 

The formula for determining the second increment of the 
correction is 



where 



[v((7,,+2A(7,,)] 



is the boundary-induced upwash 



angle corresponding to tae free-air load plus the reflection- 
plane load increment. The quantity ACi therefore corre- 
sponds to this reflection-plane load and must be divided by 
2A0,, . 



' in order to be based on unit tunnel rolling-moment 



0.C 

coefficient. 

The final corrected rolling-moment coeffident is 



1 + 



or 



0> 



1+ 



2AA, 



(20) 



2A(7,, 



where 1-| — 75—^ is obtained from figure 9 and AOjinay be 
obtained either from equation (19) or as follows: 



AC, 



(21) 



where 6[cciAlh{Cli^+2AOi^)Y is the increment of boundary- 
induced load obtained from the span-load calculations for 
the reflection-plane load Ot^-\-2AOtf. 
The effect of streamline curvatiffe was not included in 

740023 — 18 24 



equation (20) because calculations for several models showed 
that the effect was small enough to be neglected. 

The rolling-moment coefficients as computed from equa- 
tion (20)' are of the correct magnitude but apply to a wing 
angle of attack slightly different from the geometric angl^ of 
attack corrected for the symmetrical-load boundary- 
induced upwash angle determined for zero rolling moment. 
The effective change in wing angle of attack results from the 
aileron boundary-induced upwash angle and the reflection- 
plane induced angle. The fact that the corrected rolling- 
moment coefficients really apply to a slightly different wing 
angle of attack is of importance only near the stall or for 
aileron arrangements particularly sensitive to angle-of-attack 
changes and is usually neglected. The angle. -ujhange is 
usually small, less than J4°. 




for reflection- 
plane models 

in 7-by 10-foot 
closed rectangular 
wind iunnala 

a4. 



.7 

— 1.0 



A .6 
Tcpcr ratio. 



10 



12 



FiouBB 10.— Values of Uie aerodyoamlo-indnotlon factor J tued to calculate the Toning- 
moment ooneotion (or asjnnmetilcally loaded leflectlon-plane models for ttnee diflerent 
Taciatlons of bonndaiy-lnduced ap^rash angle. These npwash-angle variations corre- 
spond closely to those (or alleron^pan ratios of 0.4, 0.7, and 1.0 (or models In 7- by lO-foot 
closed reetangnlar wind tunnels. 

Induced yawing moment. — The correction to the induced 
yawing moment results from the interaction of the several 
components of boundary-induced upwash velocity and the 
several components of load as well as from the reflection im- 
age and the absence of the reflection wing. The calculation 
procedure will be to determine separately each component 
of the correction and then to sum up the various components 
as follows: 

Aa,=(AC,0i+(Aa.,).+(Aa,)3+(AC.,).+(Aa.,)5 (22) 

where the various components are defined in the symbols. 

Values of the correction due to the reflection plane 
(A(7b,)i were calculated from the influence lines of reference 5 
for a series of aileron-span ratios for the plane wing and for 
flap-span ratios of 0.4 and 0.7. TMs part of the yawrng- 
moment correction is due solely to the reflection plane and 
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does not depend upon the tunnel-model dimensions. The 
effect of flap span proved to be negligible; so the values of 
the correction as presented in figure 11 in the form of curves 
of (.AOn^JOiJJi. against aaleron-span ratio for* ailerons 
extending inward from the tips are therefore for values of 
the lift coefficient equal to the measured lift coeffident. 
The other components are determiaed from equation (13) 
where the product wT is defined by the subscripts 2, 3, 4, 
and 5. (See STjMBOLS.) All components of boundary- 
induced upwash velocity and load have already been calcu- 
lated in the form of parameters that are easily converted to 
the product wT. 

Hinge moments. — The jetr-bovrndary corrections to the 
hinge moments of lateral-control devices are usually small. 
The correction due to the elliptical streamline-curvature load 
for tiie symmetrical loading condition has already been 
presented. Another small correction exists because the load 
due to aileron deflection is greater in the wind tunnel than 
in free air. Although the load due to aileron deflection may 
be as much as 15 or 20 percent greater ia the wind tTumel, the 
correction to the hinge monients is very small because the 
greater part of the boimdary-induced load is of the additional 
type (similar to that produced by an angle-of-attack change), 



which has only a small load over the aileron portion of the 
wing. The correction is calculated from the average diffe]> 
ence between the span loading due to the lateral-control 
device in the wind tmmel and the loading in free air due to 
a given aileron deflection. Because the correction is small, 
the inerement of load at the aileron center section 

will usually be sufiB.cient for the calculation. The correction to 
the aileron hinge-moment coefficient is then assumed to equal 
the correction to the section hinge-moment coefficient of the 
aileron center section (or of some other typical section) 

where &Ca/£)c„ may be determined experimentally or from 
section data such as reference 9 for plain-flap lateral-control 
devices. The hinge-moment correction resulting from the 
aileron-load streamline curvature is usually small enough to 
be neglected. 
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Ai/ei~on-apan rafio, . 

Tmna IL— lociemcnt of yawing-moment ccmectlQn due solely to the leflectton plane tor aflennu extending Inward from tlie tips of tlie wings of reference 6. 
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NUMERICAL VALUES OF CORRECTIONS FOR MODELS IN 
7- BY 10-FOOT CLOSED WIND TUNNELS 

Some nuinerical values of the more important corrections 
were computed for various reflection-plane models mounted 
in 7- by 10-foot closed wiad tunnels. Cross plots were made 
to deteimine the variation of the corrections with each of the 
model parameters — A and X, for example — and the residts 
are presented as graphs and empirical equations that may 
easily be used to estimate the values of the corrections for 
almost any model. The computations were made by using 
the load curves and the chord distributions of reference 5 
for constant-percentage-chord flaps extending outward from 
the plane of symmetry and for constant-percentage-chord 
ailerons extending inward from the wing tip. The values of 
the corrections presented should be suflBiciently accurate for 
models that deviate slightly from these conditions. The 
corrections are given in terms of the measured lift and the 
measured rolling-moment coefficients. 

The corrections to the angle of attack, the lift, the drag, 
and the yawing-moment coefficients for symmetrical loading 
conditions are presented in figure 12 for wings having X=0.5 
and without partial-span flaps. The corresponding correc- 
tions for models of any taper ratio and flap-span ratio may 
bo dotormined by substituting the values obtained from 
figure 12 in the following enipirical equations: 

Aa=(^)^ , [l-0.063(X-0.5)](Ci;„-l-0.90(7x,p 

" /AC \ 

AC7^... =(-^)^^ [l-0.116(X-0.5)] 

AOd,=(^^^ [1 -0.080 (X-0.6)] 

{c./-|-2[l-l-0.20(l-^)]c.„C7.,+[n-0.65(l-||)](7V 
AC,,=(^)^^ [H-0.130(X-0.5)] 

The corrections to the rolling- and the yawing-moment 
coefficients for asymmetrical loading conditions may be 
determined from figure 13. The principal curves of this 

figure are drawn for models having -^=7 feet and X=0.50. 

Supplementary curves on the same figure give additional 
correction increments that account for the effect of other 
spans and other taper ratios. .The total corrections are 
obtained by adding the additional increments to the correc- 
tions obtained from the principal curves, as follows: 



-■.=[(m,.-(a).J-'-' 

The streamline-curvature corrections to the hinge-moment 
coefficients of plain or balanced (overhang type of balance) 
flaps and ailerons were computed.' The balance chord is 
assumed to be a constant percentage of the flap or aileron 
chord, and the flaps and ailerons are also of constant- 
percentage types. Figure 14 gives values of the parameter 

ACkA^(^^ /Cifi^F for various flap- and aileron-span ratios 
and various taper ratios. The fa ctor F is taken from figure 8. 

It should be remembered that the streamliae-curvatiu^e 
correction to the hinge moments is a function of the cube 
of the chord and comparatively minor variations in plan 
form, such as tip shape, thus may change the correction by 
10 or 15 percent for a given wing span and a given aspect 
ratio. The value of the correction that is determined for 
a wing with linear taper in the following fllustrative^example 
indicates the possible change in the hinge-moment correction 
with plan-form details. The correction is xisually fairly 
small, however, so changes of 10 to 15 percent in its magni- 
tude are not too critical. 

The correction to the aileron hinge-moment coefficient 
caused by asymmetrical loading may be expressed roughly as 

where the -value of (d in degrees) is determined ex- 
perimentally (or estimated) for a given model. This cor- 
rection is very small for models having ordinary proportions. 

ILLUSTRATIVE EXAMPLE 

The method of computing the jet-boxmdary corrections for 
a reflection-plane model will be illustrated in detail by an 
example. The method applies to any rectangular timnel 
but, because the computations of the boimdary-induced 
upwash velocities (figs. 4 and 5) have been performed only for 
a 7- by 10-foot closed tunnel, the example is for a model 
ia a tunnel of these dimensions. 

In practice, however, the principal corrections for models in 
7- by 10-foot tunnels can be more easily obtained from the 
graphs and empirical equations just presented. In practice, 
also, it will seldom be worth while to compute aU the cor- 
rections, such as those to the span load distribution; however, 
for completeness, all the corrections will be computed in the 
example. It might be noted that it is often convenient to use 
a single average correction for wing-flap combinations rather 
than to break the correction into two parts. The accuracgr 
required for the corrections will determine the number of 
these simplifications that may be used. 
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Aspect ratio, A 



FiooBE 13.— Charts for dotormlnlng tbo coirooUons to tho rolling- and yawlng-momont coefll- 
otonts of asymmetrically loaded ivflc«ttoii-plano models mounted in 7- by 10-foot domd 
rectangular wind tunnels. 
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FtaoBK M.— Obarta for determining tho correotlomr, duo to streamline curvature, to tho flap 
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FioOBE U^FIan fonn of the sembpan 'Ring modd used lor the UlDstiBtlve example. 
X=6; X-aS ^-(U; ^-OJL 

The constants for the assumed model (fig. 15) that are 
reqxiired for the computations are as follows: 

Aspect ratio, A 6 

Taper ratio, X a 5 

Semispan, 6/2, feet 7. 0 

Wing^area, S, square feet 32. 67 

Flap area, Sr, square feet ." L 94 

Flap-span ratio, ^ 0. 5 

Flap-chord ratio, c//c l_ 0. 2 

Flap mean chord, c/, foot 0. 65 

Aileron area, jSo, square feet L 19 

Location, inboard aileron tip, ^ 0. 50 

Location, outboard aileron tip, ^ 0. 97 

Aileron-chord ratio, cjc 0. 2 

Aileron mean chord, 5a, foot 0. 40 

The wing has rounded tips and is equipped with plain, 
unbalanced, sealed, constant-percentage-chord " flaps and 
ailerons. The model is mounted on the center line of a 
7- by 10-foot wind tunnel. 

Because the influence lines of reference 5 are used in the 
calculations, the slope of the section lift curves used for 
reference 5 will be used in this example. The value of the 
slope of the section lift curve Oq is 5.67 per radian and of the 
slope of the finite-span lift ciu*ve Oi is 4.38 per radian. 

SYAIMISTRICAL LOADING CONDmONS 

Computations for the symmetrical loading conditions 
may be made in the following steps: 

1. Values of i'w/r)„ for the horseshoe vortex representing 
the wing are obtained from table II by assuming an effective 
vortex semispan of 6.0 feet. The recommended value of 

j/i would be 0.88 or 6.16 feet, but the nearest H-foot value 

is selected in order to use the numerical values of the upwash' 
velocities from the table without interpolation. The upwash 
angle at each station for imit lift coefficient is obtained from 

equation (1) where the single load increment ^(^^^ 

is equal to — and therefore 
Vi 

{tcX=(3X LiMxe} (3X 



The upwash angle is plotted in figure 16. 
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FlOUSB 16.— Variation of the bonndary-Indnccd upwash angle, duo to the wing and du$ to the 
flap, along the lemispan o[ the model osed for the illustrat hre example. 



2. Values of (w/r) f for the horseshoe vortex representing 
the flaps are obtained from table II for a value of j/i of 4.6 
feet, which corresponds approximately to the recommended 

ratio g^=1.29. The upwash angle is 

(w>(f)iOT-6]="-(T), 

and the numerical values axe plotted in figure 16. 

3. The free-air span loading of this wing, as obtained from 
reference 6, is plotted in figure 17. The increment of load 
due to the jet boundaries is calculated from the upwash- 
angle values determined in step 1 and the influence lines of 
reference 5 and is plotted in figure 18. The area under the 
curve of figure 18 is 0.583, which corresponds to a (AOlIOz,)^ 
of 0.583/7 or. 0.0833. The increments of jet-boundary- 
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FiQUKE IZ.— Fiee^air and tnnnel span load distribatlons for the wing used for the iUustrativo 

example. 
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FIOUBB 18,— Variation of tbo boimdary-liidaoed Increment o( load, doe to the wing and due 
to tbe flap, along tbe lembpan of the model used for the niustratire example. 

induced load are multiplied by 1/(1—0.0833) or 1.091 and 
are added to the free-air load and plotted in figure 17. The 
tunnel loading corresponding to imit lift coefficient is 
obtained by multiplying the ordinates of the total-load 

7 0 

curve by the area ratio 7 q-^q 5^3x1 qqi ' ^ comparison 
of the final tunnel load curve with the free-air load curve 
and the span-load data of reference 6 shows that the tunnel 
loading corresponds to a tape ratior of about 0.45. 

4. The same procedure applied to the span-load curve 
for the flaps, as obtained from reference 5, is illustrated in 
figures 18 and 19. 



5. The values of 



da; 



required to compute the cor- 



rections due to streamline curvature axe token isom table I. 
The summation product 



dx 



1.6 



1.2 




5 



rree-air , . _ .i 

Tunnel \Q.*AQ./k-^)\ 
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FioUBB 10.— Flre»alr and tunnel span load distributions for the flap of the model used for the 
fliustratlva example. 



for the wing corresponds to a one-step approximation to' the 
load curve as shown in figure 20. The value of 

A for yi=6 is 1.167. The results of the calculation 

of the desired product 

are given in figure 21. 
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W Whig, 
(b) Hap. 

TroQBE 20.— One^p approximation to the wing load onrve and two4tep approximation to 
the flap load curve used for oompntatlon of coneotions due to streamline currature for the 
illustrative example. 



6. The product 



6 L cte 



for the flap is found 



_l t 



similarly. The two-step approximation to the load curve is 
shown in figure 20 and the numerical values of the final 
product are given iu figure 21. 
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or, from step 7b, 
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Dletanco from plane of symmetry, y, ft 
FiQUBB 21.-Varia«on of M[a£2£^-j^^ ^i^.^JS^-^^^ atoDg the »mW 
' of the model used for tbe IIlaztTative example. 



7a. The correction to the angle of attack due to tiie effect 
of the jet boTindaries on the mng alone is the sum of the 
corrections obtained from equations (3) and (6). The value 
of (AC!t/Ci)o in equation (3) is 0.0833 as found in step 3. 
Thus, from equation (3), 

. _ 0.0833X57.3 ^ .q.^ 



The factor 0.50- 



of equation (5) may be assumed to be 



equal to 0.25 for the wing alone and 



X 



m2A 



'(m) 



(? (iy=0.248 



from a mechanical integration of the curve in figure 21. 
(The moment of the area about the plane of symmetry is 
found at this time to be 0.580 and will be used in step 14.) 
Thus 



57.3 X0.25Q; „^ 
14 



'0.248=0.254a;„ 



7b. The alternate method of obtaining the angle-of-attack 
correction is given in equation (4). Mgure 22 gives values of 

(^^^"^ c for the wing. The area under the curve is equal to 

0.307, which gives a correction equal to 

57.3a,„X2X0.307 
Aatt,= 32 67 =1.078Oi;„ 

The value of (Aa:,.e Jif wiU,. of course, be the same as found 
in step 7a. 

8. The total angle-of-attack correction due to the wing 
alone is, from step 7a, 

Aa„= 1.090C7i„-{-0.254a;„= l'.344a;„ 



A«„= 1 .078Ci;„-l- 0.264(7£„= 1 .332Ct„ 

9. The correction to the angle of attack due to the ilap 
loading is found, from equation (3), to be 

Aay=1.200Ci, 

or, from equation (4), 

Aa/=1.190(7£, 

The effect of streamline curvature is not considered because 
the factor 0.50— —-■ft'O for the flapped wing. 

10. The total correction to the angle of attack is the sum 
of the wing and flap corrections, from equations (3) and (6) 

Aa,o^,= Aa„-t- Aa^= 1.34ai„-|- 1 .20(7i, 

or, from equations (4) and (6) 

Aa,.„,=1.33ai„+1.19Ct/ 

11. The correction to the lift of the wing is given in 
equation (7) as 



-0.0194(7, 



12. The correction to the lift of the flap is computed from 
equation (6) as 

, ^ 4.38(7i, 

i^Or.,,,), 4X14-0-267 0.0209(7., 

13. The total lift correction is 

(ACi,..c.).»r.,=-0.0194a,„-0.0209(7^, 
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FiooBE 23.— Vartatlon of c along the semlspan of the model osod for the niustrallvo 

example. 
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FinOBG 23.— Variation of tbe parameten used to compate the indnced-drag and yaving- 
momont oorrcotloDS, Ua symmetrical loading, along tlie semlspen o{ tlie model used Ibr the 
tUostratlTO example. 

14. The moments of area determined in steps 7 and 12 
will be used in the determination of the spanwise center of 
pressure of the lift correction load from equation (9), <?hich 
gives a value of y,.c.=2.36 feet for the mng correction and 
yf.e.=2.16 feet for the flap correction. These values axe 
required for the computation of the correction to the 
pitching-moment coefficient. 

15. The correction to the measured pitching moment due 
to the wing alone is obtained from equation (8) where the 
chord c at the lateral center of pressure determined in step 
14 is 2.64 feet. Because the pitching-moment coefficients 
for this model are given about the 0.25c line and in terms 

2 64 S 

of the mean chord, arc„=-^ and Cu,=-j=2.333. The correc- 
tion as obtained from equation (8) is 

..^ . 2.64 5.670,, 
(.^t;m...J„- 4 2.333 X4X14"-^*^ 

=0.007iai„ 

16. The correction to the flap pitching moment is found 
by the same method to be 

(^C'„....),=0.0077Cii, 

17. The total pitching-moment-coeffident correction is 

Aa„.. , =0.007iai„+0.0077C!i;^ 

* 18. The induced-drag-coefficient correction is determined 
from equations (11) and (12) and the integrals of the various 
products of boundary-induced upwash and load parameters 
are obtained in steps 1 to 4 and are plotted in figure 23. 



Thus, 

A(7b,=3| [0-132CD„^-l-0.169(7i/-l- (0.145+0.144)Ci„<7i/] 
=0:0189Ci,„^-|-0.024lCi,/+0.0413C7z„C7£, 

19. The induced-yawing-moment-coefficient correction, 
obtained from equations (12) and (13) and the moment 
integrals of figure 23, is 

AO.,' [0.340(7i„^-|-0.306(7i/-f- (0.286+0.330)(7i„<7i^] 

=-0.0017C7iJ'-0.0016ai/-0.003lCi;„Ci;^ 

20. .The correction to the flap and aileron hinge-moment 
coefficients due to streamline curvature is given in equation 
(17). The value of F from figure 8 is 0.036. Mechanical 
integration of tiie curves of figure 24 gives for tiie corrections 
to the flap hiage-moment coefficient 

. ^ 5.67X14X0.036 , ^ , 

^^V=8^X6X1.94X0.55 (0.602(7x„+0.602C!.,) 



or 



ACA^=0.0089a;„+0.0107a;^ 



The value of the correction detemuned from figure 14 is 
somewhat smaller (A(7»^=0.0083(7i), because the chords near 
the root sections are smaller for the wings of reference 5 
than for wings with linear taper. The correction to the 
aileron hinge-moment coefficient is 

^ _ 5.67X14X0.036 , ^ 4.0 ^ 
^*^*«~87rX6X1.19X0.40^"-^^^^^''+"'^^^^^^^ 

A(7s„=0,005ia;„+0.0046a;.-. 
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FIOTOB M.-Variatlon of^ ^J^£|^ ] .C and ^ j- SWVCt) j^^ ^long the semlspan of 
tbe model nied for tbe illnztratlve ernmple. 
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ASYMMETBICAL LOADING CONDITIONS 

21. The free-aif span load distribution due to the deflec- 
tion of one aileron is obtained from reference 5 in terms of 
cc,/c,a and is plotted in figure 25(a). The loading in terms 
of the parameter cCiA/bOi^ is obtained by so adjusting the 
ordinates of the curve of figure 25(a) that the moment of 
area about the plane of symmetry is equal to 4.0(6/2)*. The 
conversion may be made graphically or by means of equation 

(18) where ^^=1.37 from figure 6 and, from figure 16 of 

reference 5, the value of Oijk is the difference between the 

values at ^=0-^7 and at ^=0-5; that is, 

-^=0.73—0.25=0.48 



and 

Therefore 



^=0.5^=0.24 
a k 



bOi. 0.24 C,a ' cot. 



The new curve of free-air load is plotted in figure 25(b). 
The reflection load is added to the free-air load to give the 
reflection-loading curve of figure 25(b). The jet-boundary- 
induced upwash angle is obtained from equation (2) for the 
three-step approximation to the load curve (as indicated in 

fig. 25(b)). The nximerical values of [^ ^((j^ ^2 AO, ) j 

plotted in figure 26. The corresponding increment of load 
at each spanwise station is obtained from the influence lines 
of reference 5 and is presented in figure 27. The increment 
presented in figure 27 must be divided by 



1 + 



2Ao: 



whore A(7{ is obtained from equation (21) and the moment of 
the curve of figure 27 as 



and 1-f 



2Ag,, 



AC,=iY^=0.0547 

is equal to 1.084 from figure 9. The corrected 



increments are added to the reflection-load curves of figure 
25 (b) to give the tunnel-load curve (fig. 25 (c)) at the same 
aileron angle as the curve for the free-air load. The ordinates 
are again adjusted to give a' moment of 4.0(6/2)*, which 
corresponds to unit rolling-moment coefficient. The result- 
ing corrected tunnel-load curve is presented in figure 25 (c). 

22a. The corrected rolling-moment coefficient is obtained 
from equation (20) where 



2A(?,, 



=1.084 



and, from step 21, 



A<7j=0.0547 



0.0547 \ 

" 1.084 ; 



=0.876(7, 



1.084 

22b. The alternate method of determining the corrected 
rolling-moment coefficient is to use equation (19) to calculate 

AC/. The product -^2AG )J " ^ plotted in figure 28 

and the area moment is found by mechanical integration to 
be 5.848. Figure 10 gives J=1.93. Therefore, 



AC,= 



6*X5.67 



14^X(6-fl.93) 



X5.848=0.0547 



which is the same, in this case, as the value calculated in 
step 21. The agreement usually will be close but not 
necessarily exact. 

23. The correction to the yawing-moment coefficient duo 
to aileron deflection is obtained from equation (22). Tor 
thia model the value of (A(7n,)i, due to the absent wing, is 
foimd from figure 11 to be — 0.0104(7/£Ci. Step 22a gave 
O}=0.876Ot. Thus: (AO,^, =—0.00900,(71,. By consider- 
ing the afleron-wing combination and using (w/VOi,)a from 

figure 16, ^ y^^^ +2AOi ) J from figure 26, tunnel (cc,.4/6Cx,)» 

from figure 17, and tunnel cctAlbGig from figure 25(c), 
equation (13) gives 



(AC,.),= - 



and 



dy 



The aileron-flap combinations are determined similarly. 

The products of the various upwash and loading parame- 
ters are plotted in figure 29 and mechanical moment integra- 
tions give 

2.442 



(AO.,), 



(AO..), 



2.925 



J2 Ol^Ci 

1.680 
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(a) Free air loading in temu of Data from relaonce S. 



(b) FtokIt loading and reflection loading in tenna of lre»air pins reflection loading in tarnn of j^; ,^^^^ +2 1<?^ 3 * 
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FiaOBE 25.— Devdopment of the tminel span load distrlbntlon due to tbe deflection of tbe aileron of the model used for Che llliutrative example. 
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From reference 9 

|^*=--0.265^ 
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Fin DBS 28.— Variation of tUe boimdax:r-Sndaeed upnash aii(^ doe to aDeron deflection, 
alons the semtepan of the model used for tbe IDnstratfve """t^" 



2A(7 

When it is known that 6=14 feet and 1+ >, '^ =1.084, the 

total correction to the yawing-moment coefficient as obtained 
from equation (22) is 

A(7,<=-0.009CiC,-0.026Ci„C,-0.02lCi,^C, 

25. The second component of the aileron hinge-moment 
correction is obtained approxtmately from equation (23), 
with the use of figures 27 and 25 (b) . The total increment of 
load at the aUeron center section is 

^ [KC7,+2l>,)] +<^)r°-"^+°-^^=°'^° 
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FlQUBE 27.— Variation of the boandaiy-indnoed increment of load, doe to allenm deflection, 
along the scmlspan of the model osed for the UlnstrattTe example. 



Therefore, the correction is 



^^*'^~ 2.0X6.0 (-^'•^SS) (0.2) (0.70) =0.04C7, 
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FiQDBK 28.— Variation of +2ACi^ l ' semispan of the model used tof tho 

miistiatlTa example. ^ 

CONCLUDING ' KEMAEKS 

The method for determining the jet-boimdaiy corrections 
for reflection-plane models in rectangular -wind tunnels was 
presented in some detail in order to make tlie method as 
routine as possible. The method includes the determination 
of the tunnel span load distribution and the derivation of 
equations giving the corrections to the angle of attack,, 
the lift and drag coefficients, and the pitching-, rolling-, 
yawing-, and hinge-moment coefficients. The principal 
effects of aerodynamic induction and the curvature of 
the streamlines have been considered. 

Nimierical values of the more important corrections wore 
calculated for a series of representative models mounted 
in 7- by 10-foot closed rectangular wind tunnels. In order 
to simplify the calculation of corrections for models of 
uniisual size in 7- by 10-foot closed wind tunnels, tables 
of the nimierical values of the jet-boundary-induced upwash 
were presented. 
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FiODBE 29,— Variation of the porometers used to detennine the Indaoed-yawing-inomeiiit 
correction duo to aileron deflection along tbe <MTni«paii ol the model used Ibr the 
Ulustrotivo example. 
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TABLE I 

BOUNDARY-INDUCED UPWASH VELOCITY BEHIND LIFT- 
ING LINE DUE TO SINGLE UNIT COUNTERCLOCKWISE 
VORTEX ON TUNNEL CENTER LINE AND AT TWO DIS- 
TANCES yi FROM REFLECTION PLANE IN 7- BY 10-FOOT 
CLOSED RECTANGULAR WIND TUNNELS 



\ 

\ X 

S \^ 


0 


0.6 


1.6 


3.0 


6.0 


9.0 


atx-0 




0 
2 
i 
6 
8 


a 01325 
.01096 
.OOISS 
.00146 
-.00037 


0.01633 
.01277 
.00700 
.00194 
-.00018 


0.01936 
.01616 
.00900 
.00281 
.00024 


a02397 
.02003 
.01146 


0.02S06 
.02361 
.01363 


a02S61 
.023»t 
.01380 


0.00412 
.00362 
.00222 
.00096 
.00036 


iri-8 


0 
2 
4 

e 

8 


a 01706 
.01606 
.01286 
.00827 
.00661 


0.01997 
.01880 
.01621 
.00994 
.00666 


a02641 
.02401 
.01967 
.01296 
.00882 


0.03193 
.03018 
.02666 


0.03787 
.03613 
.02923 


a03821 
.03618 
.02962 


0.0OSSS 

.co«o 

.00472 
.00334 
.00230 
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TABLE n 

BOUNDARY-INDUCED UPWASH VELOCITY AT THE LIPTING LINE DUE TO A SINGLE UNIT COUNTERCLOCKWISE VORTEX 
AT VARIOUS DISTANCES yi FROM THE REFLECTION PLANE IN 7- BY lO-FQOT CLOSED RECTANGULAR WIND TUNNELS 





as 


1.0 


L6 


2.0 


2.6 


3.0 


3.S 


4.0 


4.6 


6.0 


6.6 


e.0 


6.fi 


7.0 


7.8 


8.0 


S.8 


9.0 


S.6 


V (ft) \^ 









































c 

.6 

LO 
L6 
2.0 
2.6 
3.0 
3.8 
4.0 
4.6 
6.0 
6.6 
6.0 
0.6 
7.0 
7.6 
8.0 
&S 
9.0 
9.6 

lao 



0 
.6 
1.0 
LS 
2.0 
2.6 
3.0 
3.6 
4.0 
4.5 
8.0 
6.8 
6.0 
6.8 
7.0 
7.5 
8.0 
8.5 

ao 

9.S 

mo 



Vortex on tannd center line 



0.00267 
.00202 
.00249 
.00229 
.00202 
.00173 
.00142 
.OOUl 
.00034 
.00081 
.00038 
.00018 
.00006 
-.00004 
-.00014 
-.00018 
-.00021 
—.00022 
-.00018 
-.00013 
-.00006 



a00287 
.00282 
.00260 
.00241 
.00210 
.00176 
.00141 
.00103 
.00078 
.00053 
.00031 
.00013 
-.00001 
-.00011 
-.00018 
-.00023 
-.00025 
-.00025 
-.00022 
-.00017 
-.OCOlO 



0.00523 
.00514 
.00490 
.00452 
.00402 
.00344 
.00284 
.00228 
.00172 
.00122 
.00079 
.00044 
.00014 
-.00008 
-.00022 
-.00035 
-.00040 
-.00039 
-.00035 
-.00024 
-.00011 



a00764 
.00761 

.oono 

.00663 
.00594 
.00513 
.00428 
.00345 
.00264 
.00190 
.00128 
.00076 
.00030 
-.00004 
-.00029 
-.00044 
-.00053 
-.00063 
-.00046 
-.00033 
-.00010 



aoooso 

.00966 
.00924 
.00858 
.00774 
.00678 
.00574 
.00466 
.00383 
.00270 
.00188 
.00114 
.00057 
.00009 
-.00026 
-.00047 
-.00057 
-.00059 
-.00051 
-.00031 
-.00004 



0.0U67 
.01153 
.01108 
.01035 
.00942 
.00836 

.oono 

.00592 
.00472 
.00359 
.00256 
.00168 
.00093 
.00036 
-.00009 
-.00039 
-.00053 
-.00055 
-.00046 
-.00022 
.00018 



a 01325 
.01310 
.01284 
.01192 
.01096 



.00853 
.00722 
.00588 
.00458 
.00341 
.00236 
.00146 
.00075 
.00022 
-.00015 
-.00037 
-.00039 
-.00028 
.00004 
.00051 



a 01462 
.01437 
.013» 
.01325 
.01230 
.01114 



.00349 
.00708 
.00570 
.00137 
.00319 
.00217 
.00133 
.00069 
.00024 
-.00001 
-.00008 

.omio 

.00047 
.00105 



a 01548 
.01636 
.01498 
.01432 
.01343 
.01238 
.01110 
.00972 
.00831 
.00687 
.00548 
.0(M19 
.00306 
.00211 
.00135 
.00083 
.00053 
.00048 
.00065 
.00111 
.00188 



a 01620 
01609 
01674 
01616 
01438 
01339 
01222 
01092 
00951 



.006^ 
.00535 
.00413 
.00308 
.00225 
.00164 
.00132 
.00128 
.W149 
.00201 
.00298 



a 01670 
01658 
01627 
01680 
01612 
01424 
01321 
01201 
01070 
00933 
00796 
00683 
00537 
00427 
00337 
00274 
00237 
00233 
00287 
00338 
00453 



Ol 01698 
01690 
01684 
01623 
01588 
01494 
01403 
01299 
01183 
01057 
00927 
00798 
00677 
00568 
00478 
00410 
00375 
00378 
0O423 
00514 
00061 



0.01TO8 
01702 
01684 
01850 
01606 
01645 
01472 
01385 
01288 
0U77 
01059 
00941 
00827 
0O728 
00639 
00577 
00551 
00584 
00825 
00745 



a 01708 
01702 
01688 
01666 
01629 
01683 
01627 
01459 
01379 
01288 
01191 
01088 
00990 
00900 
00827 
00780 
00768 
00798 
00887 
01049 
01304 



0.01698 
.01694 
.01684 
.01667 
.01644 
.01611 
.01670 
.01521 
.01461 
.01393 
.01317 
.01340 
. 01161 
.01091 
.01041 
.01018 
.01027 
.01089 
.01222 
.01446 
.01795 



0.01630 
01680 
01673 
01062 
01649 
01631 
01605 
01572 
01636 
01490 
01442 
01390 
01341 
01302 
01280 
01288 
01339 
01451 
01848 
01968 
02471 



a 01662 
01659 
01658 
01665 
01649 
01843 
01833 
01619 
01601 
01584 
01663 
01643 
01631 
01630 
01549 
01603 
01713 
01898 
02197 
02872 
03441 



0.01838 
01040 
01641 
01645 
01649 
01651 
01657 
01662 
01668 
01074 
01686 
01704 
01732 
01778 
01853 
01873 
02102 
0«6S 
02902 
03070 
04973 



0.01018 
01620 
01027 
01636 
01047 
01063 
01680 
01700 
01736 
01709 
01816 
01874 
01951 
02066 
02202 
02113 
03710 
03188 
03931 
05221 
07854 



Vort«i 0.6 foot above or below center Uno 



0.00563 


0.00818 


a 01044 


a 01237 


a 01396 


a 01520 


0.01612 


0101978 


0.01718 


a 01739 


0. 01744 


0.01737 


a 01721 


0.01699 


0.01675 


a 01649 


a 01026 


.00552 


.00803 


.01027 


.01220 


.01378 


.01504 


.01598 


.01665 


.01703 


.01731 


.01738 


.01732 


.01718 


.01698 


.01674 


.01650 


.01628 


.00522 


.00762 
.00698 


.00979 


.01168 


.01328 


.01457 


.01656 


.01629 


.01878 


.01707 


.01719 


.oin9 


.01709 


.01693 


.01674 


.01063 


.01633 


.00175 


.00903 


.01087 


.01247 


.01380 


.01488 


.01570 


.01628 


.01668 


.01688 


.01696 


.01691 


.01685 


.01671 


.01050 


.01643 


.00417 


.00617 


.00806 
.00696 


.00982 


.0U40 


.01278 


.013»4 


.01486 


.01558 


.01600 


.01643 


.01663 


.01672 


.01672 


.01668 


.01061 


.01088 


.00351 


.00525 


.00859 


.01013 


.01153 


.oi2n 


.01382 


.01468 


.01636 


.01584 


.01619 


.01642 


.01655 


.01663 


.01660 


.01072 


.002S4 


.00430 


.00578 


.00726 


.00872 


.01012 


.01142 


.01258 


.01350 


.01443 


.015U 


.01562 


.01602 


.01631 


.01653 


.01073 


.01093 


.00220 


.00337 


.00461 


.00591 


.00725 


.00880 


.00993 


.01U8 


.01233 


.01336 


.01421 


.01493 


.01652 


.01600 


.01642 


.01679 


.01718 


.00161 


.00251 


.00350 


.00460 


.005<9 


.00708 


.00837 


.00968 


.01094 


.01213 


.01317 


.01411 


.01493 


.01663 


.01626 


.01687 


.01747 


.00110 


.00174 


.00250 


.00339 


.00441 


.00556 


.00682 


.00813 


.00916 


.01077 


.01201 


.01316 


.01422 


.01518 


.01008 


.016» 


.01784 


.00066 


.00109 


.00163 


.00231 


.00316 


.00116 


.00532 


.00660 


.007£6 


.00938 


.01076 


.01212 


.01343 


.01466 


.01588 


.01706 


.01829 


.00030 


.00054 


.00090 


.00140 


.00206 


-.00292 


.00394 


.00515 


.00649 


.00794 


.00947 


.01102 


.01267 


.01410 


.01661 


.01720 


.01880 


.00002 


.00011 


.00031 


.00065 


.ooue 


.00185 


.00274 


.00384 


.00513 


.00660 


.00820 


.00991 


.01170 


.01354 


.01544 


.01744 


.01901 


-.00020 


-.00022 


-.00014 


.00007 


.00044 


.00098 


.00174 


.00273 


.00395 


.00539 


.00705 


.00888 


.01039 


.01304 


.01635 


.01788 


.02001 


-.00038 


-.00046 


-.00046 


—.00038 


—.00010 


.00033 


.00007 


.00186 


.00299 


.00440 


.00607 


.00803 


.01023 


.01389 


.0I5M 


.01852 


.02203 


-.00044 


-.00059 


-.00066 


-.00063 


-.00016 


-.00012 


.00014 


.00123 


.00230 


.00367 


.00537 


.00742 


.00983 


.01283 


.01686 


.01902 


.02106 


-.00047 


-.00065 


-.00076 


-.00077 


-.OOOM 


-.00035 


.00014 


.00089 


.00191 


.00327 


.00501 


.00718 


.00982 


.01300 


.01681 


.02139 


.02704 


-.00046 


-.00065 


-.00076 


-.00078 


-.00066 


-.00038 


.00010 


.00082 


.00188 


.00325 


.00508 


.00742 


.01034 


.01400 


.01853 


.02423 


.03102 


-.00042 


-.00057 


-.00068 


-.00065 


-.00051 


-.00021 


.00030 


.00107 


.00217 


.00387 


.00560 


.00825 


.01169 


.01538 


.02142 


.02870 


.03899 


—.00032 


-.00043 


-.00046 


-.00040 


-.00020 


.00016 


.00076 


.OOIM 


.00288 


.<mB7 


.00684 


.00983 


.01378 


.01001 


.02611 


.03618 


.05184 


-.00019 


-.00021 


-.00017 


-.00001 


.00028 


.00077 


.00151 


.00257 


.00104 


.o&m 


.00875 


.01237 


.01725 


.02401 


.03377 


.01919 


.07814 



